Localization of cold atoms in state-dependent optical lattices via a Rabi pulse 
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We propose a novel realization of Anderson localization in non-equilibrium states of ultracold 
atoms in an optical lattice. A Rabi pulse transfers part of the population to a different internal 
state with infinite effective mass. These frozen atoms create a quantum superposition of different 
disorder potentials, localizing the mobile atoms. For weakly-interacting mobile atoms, Anderson 
localization is obtained. The localization length increases with increasing disorder and decreasing 
interaction strength, contrary to the expectation for equilibrium localization. 
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The control and manipulation of ultracold atoms in 
magnetic and optical traps has made enormous pro- 
gresses towards the realization of fundamental condensed 
matter phenomena ^, [2] ■ Among them, the phenomenon 
of quantum localization of matter waves induced by a 
random potential, observed in a variety of condensed- 
matter setups 3:], takes a new dimension in the context 
of cold atoms. In fact the capability of fine tuning the 
disorder potential as well as the interaction among atoms 
within the same sample leads to the possibility of in- 
vestigating the Anderson localization transition with un- 
precedented control, and to monitor the complex effect 
of interparticle interactions ^4] . Substantial progress has 
been made recently in this direction with optical disorder 
potentials, leading to the observation of Anderson local- 
ization for weakly-interacting bosons [5|-|7[ and to exper- 
iments aiming at the realization of a strongly interacting 
Anderson insulator (Bose glass) with bosons in optical 
lattices ^E^. 

An alternative route to Anderson localization of cold 
atoms consists in creating a disorder potential by pinning 
a random distribution of an atomic species b. A second, 
mobile atomic species a experiences a disorder potential 
due to interspecies collisions. This localizes species a 
14| . This approach has the unique advantage that 
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the disorder potential can be prepared in a quantum su- 
perposition of different realizations, such that the time 
evolution of the mobile a-atoms samples all realizations 
of the potential at once, and measurements on the mobile 
particles deliver disorder-averaged observables 12, llJ] ■ 

In this paper, we theoretically explore a different 
scheme leading to localization in atomic mixtures, based 
on "dynamical doping" of frozen 6-atoms into a cloud 
of a-atoms. a and 6-atoms represent two internal states 
of the same atomic species. Like Refs. J12h14| . we use 
a state-dependent optical lattice to freeze the 6-atoms 
and keep the a-atoms mobile. All N atoms are ini- 
tially prepared in state a), and occupy the motional 
ground state. Subsequently, a short radio-frequency Rabi 
pulse is applied, creating an internal-state superposition 
cos(a/2)|a) -I- sin(a/2)|6), where a is the pulse area. We 
show that the frozen 6 component creates a quantum 
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FIG. 1: Sketch of atom localization in state-dependent op- 
tical lattices. The a-atoms, initially prepared in the ground 
state of the lattice, are converted to a quantum superposition 
cos(Q/2)|a) -I- sin(a/2)|fo) by a Rabi pulse, where the 6-atoms 
are frozen by the lattice. The state of the system can be 
viewed as a-atoms (red) moving in a quantum superposition 
of random potentials created by the frozen b-atoms (blue). 



superposition of disorder potentials (see Fig. [1}. The 
wavefunction of weakly-interacting a-atoms in a two- 
dimensional lattice is studied numerically and evolves 
into a quasi-stationary profile with exponential tails. 
While we concentrate our discussion on bosons, our 
scheme should work equally well for fermions. The idea 
of dynamical doping can also generate new perspectives 
for other disordered systems. 

The analysis of our scheme is interesting for several 
reasons: First, the fact that dynamical doping results in 
localization is not obvious, because the Rabi pulse trans- 
fers a fraction sin^(a/2) of the population at each site 
into state |6). Without measuring the 6 distribution, no 
quantum projection noise is introduced and it is not a 
priori evident why the 6-atoms should generate disor- 
der. Second, it is not obvious how decoherence in the 
superposition between states |a) and |6) will affect lo- 
calization. Third, we show that the resulting disorder 
potential has negligible spatial correlations, unlike the 
most recent realizations in cold atom experiments [5l4lOl| . 
Fourth, the experimental implementation of this scheme 
should be significantly simpler than the previously pro- 



posed schemes for disorder in atomic mixtures 12l4l4|. 
These schemes involve either the diabatic preparation 
of metastable states, resulting in a disorder realization 
that is not deterministically controlled by the prepara- 
tion scheme [13|; or the individual preparation of each 
species on disjoint optical lattices, and the subsequent 
relative motion of the two lattices to achieve coupling 
of the mobile particles to the disorder potential [12, llJ] ■ 
Fifth, due to the quantum parallelism, single-shot experi- 
ments with large atom numbers directly provide disorder 
averages. 

While our proposal is quite general, we wish to point 
out that an experimental implementation is feasible with 
existing technology in ultracold gases, in a way that 
all relevant parameters can be adjusted independently. 
One could, e.g., choose ^^Rb atoms with hyperfine states 
\a) = \F = l^mp = 1) and |&) = |F = 2,mF = 2). 
For circularly-polarized lattice light, there is a so-called 
magic wavelength near 787 nm where the a-atoms see no 
lattice potential 15| . The ratio of the tunneling ampli- 
tudes Jo/ Jb can be adjusted over a wide range by tuning 
the wavelength A near the magic wavelength. Hence, 
one can make Jf, very small, thus freezing the 6-bosons. 
A Feshbach resonance [l6| for state \a) can be used to 
adjust the interaction strength Ua. Due to the state- 
selective character of Feshbach resonances, this leaves the 
interspecies interaction strength Uab unchanged. The in- 
tensity of the lattice light controls Uab/ J a- Finally, this 
choice of states | a) and | b) is stable against spin exchange 
and it is only subject to very slow dipolar relaxation [17 1 . 

The model - Our system consists of a and &-bosons 
trapped in the lowest band of an optical lattice, described 



by Uab = T-La 
Hamiltonian 



-nb+UabT.^ 



with the Bose-Hubbard 



-Ha = -Ja^ (ala, +H.c.)+^ ^<«- 1) +V',< 

(1) 
for the a-bosons, and the analogous T-Lb for the 6-bosons. 
Here a, a^ are destruction/creation operators for the a- 



bosons, n" = a^a, Vi 



m 2 



fi—i'o) is an overall atomic 



trap of frequency w/27r, and m is the atomic mass. We 
consider Uab > 0. 

The system initially contains only the a-bosons in the 
ground state of the Hamiltonian Ua, which can be ex- 
panded in the basis of lattice Fock states 
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subject to the constraint J^i ^i = -^■ 

Under the assumption that the Rabi pulse is fast com- 
pared to the Hamiltonian dynamics, the application of 
a pulse instantaneously rotates the bosonic operators at 
each site, a] — )• cos(a/2)a| -I- sin(a/2)6|. The resulting 



state can be expanded as a quantum superposition of 
many partitions of the N atoms over species a and b 

\^ab)t=,^ J2 ^(K}:K})IK})®IK}> (3) 



under the constraint ^^ 
be rewritten as 
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N. This state can 



\^ab)t=0 = E <{<})\{<}) ® \<^a{{n1}))t^0, (4) 

under the constraint ^^ n^ < N, where 
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is the normalized state of the a-bosons subject to the 
constraint of having a Fock state |{nj}) of the 6-bosons. 
c({n^}) is a normalization constant and can be chosen to 
be real. 

The subsequent time evolution of the system is gov- 
erned by Hab- The evolved state at time t, \'i'ab)t, can 
again be expanded in the form of Eq. ([3]). We now assume 
that tunneling of 6-bosons is negligible, Jb — 0, so that 
the &-bosons are frozen in the Fock state |{?T.i}). Hence, 
the a-boson components in Eq. ([5]) evolve decoupled 
from each other as \^a{n'l}))t = exp[—i'Hab{{n''i})t/h] 
\^a{{n^})}t=o- At time t, observables Ag, referring to 



a-bosons only, take the expectation value 12 1 
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which describes the average over parallel quantum evo- 
lutions of the a-bosons in all possible realizations of the 
static potential Vi = Uabn\ generated by the frozen b- 
bosons. These realizations follow the statistics dictated 
by the |c({n^})p coefficients, which are completely con- 
trolled by the choice of the initial state of the a-bosons. 
In particular, if this state is homogeneously spread over 
a given region of space, the distribution of the resulting 
effective potential will be dominated by configurations in 
which a number ~ N sin (Qf/2) of 6-bosons is randomly 
scattered over this region, hence giving rise to an effective 
disorder potential. 

With Jb = 0, the Hamiltonian 'Hf, obviously con- 
tributes only a phase factor to the time evolution of each 
\'^a{n\}))f This has no effect on the expectation values 
in Eq. ([S]). Hence, T-Lb can be omitted in our numerical 
calculations below. Furthermore, if there is spin decoher- 
ence, i.e., if the relative phase of states a and b fluctuates 
due to experimental imperfections, this will cause fluctu- 
ations in the phases of ah the \^a{n'l}))t in Eq. (JH). This 
also has no effect on the expectation values in Eq. ([6]). 
This insensitivity to spin decoherence is a strength of our 
scheme. 



Disorder statistics - In the following, we consider 
weakly interacting {Ua ^ Ja) Bose gases. To analyze 
the statistical properties of the disorder, we describe the 
state before the Rabi pulse with a condensate wavefunc- 
tion on a lattice pj] 
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which is a good approximation far from the superfluid-to- 
Mott-insulator transition (occurring for Ua/ Ja ~ 16 with 
n = 1). Here, -0^ — '^i/\'N is the effective single-particle 
state. After the Rabi pulse, the state can be written in 
the form of Eq. (U) with |$a{?^■})t=o = KV'*})®^^^^''^ 
and 
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The |cp coefficients describe the statistics of a random 
potential with an average number of frozen bosons per 
site {n\) — {Nb) |-0i| and an average total number of 
frozen bosons {N^) = A^sin^(a/2). The correlations in 
the disorder potential can be described by the density- 
density correlation function 
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Therefore, for large N the disorder potential at different 
sites is uncorrelated. Thus in contrast to speckle poten- 
tials, our proposal does not suffer from problems with 
the autocorrelation in the disorder potential. The disor- 
der distribution is, however, inhomogeneous because of 
the external trapping potential. 

Numerical simulations - In our simulations, we take 
A^ = 50 bosons, a two-dimensional LxL lattice with L — 
60 and ^uj'^<f = 0.004J<j which yields a density (n) < 1 
throughout the trap ((n) ~ 1 for the Ua = Q case). Here, 
d = A/2 is the lattice spacing. We determine \ijji\^ by 
numerically minimizing the Gross-Pitaevskii (GP) energy 
functional |17|] subject to the condition that the total 
number of particles be fixed to N . Next, we generate 
disorder realizations {n\} according to the |cp statistics 
of Eq. ([5]) via standard Monte Carlo Metropolis sampling 
14| . Finally, we numerically solve the time-dependent 
GP equation with T-Labi omitting Hb as discussed above. 

Fig.[5]shows the time evolution of the density profile of 
the a-bosons for Ua = 0, Uab = 7Ja, a = 7r/2. After the 
Rabi pulse at i = 0, the a-boson cloud rapidly expands 
from its initially Gaussian state, due to the repulsion 
from the dynamically doped 5-bosons. The expanding 
cloud wings are then refiected by the trap. This gener- 
ates breathing oscillations, which are quickly damped by 




FIG. 2: Time evolution of the density profile n{x) = {n'^/j) 
for the non-interacting gas, Ua = 0, after a pulse of duration 
a — 7r/2 and with interspecies interaction Uab ~ 7Ja- 



the disorder. Remarkably the resulting, quasi-stationary 
density profile features a central peak with exponentially 
decaying tails, (ui) ~ exp(— jr.; — ro|/^), clearly exhibit- 
ing Anderson localization. A similar behavior is observed 
for Ua > and for different values of Uab and a. Yet, the 
dependence of the localization phenomenon on these pa- 
rameters is nontrivial and differs significantly from what 
one expects at equilibrium. 

The fast, diabatic Rabi pulse takes the system out of 
the initial ground state into an excited non-equilibrium 
state. According to Eq. ([S]), the remaining a-bosons oc- 
cupy a superposition of different condensate states with 
variable particle number N — N},; these condensate states 
are non-equilibrium ones because they are immersed in 
the disordered background of frozen 6-bosons, and be- 
cause, in presence of finite interactions Ua > 0, the con- 
densate wavefunction is no longer minimizing the inter- 
action energy for a reduced particle number N — Ni,. To 
quantify the distance from equilibrium, we evaluate the 
energy transfer per a-boson due to the Rabi pulse 

r. _ ('Ha + C^a6 2_/i "-f"-!/ after {Ha)hofoTe /,„n 



A^cos2(a/2) 



N 



where (...)bcforo(aftor) denotes the expectation value over 
the wavefunction before (after) the pulse. Fig. [3^/b 
clearly shows that the transferred energy 6e increases 
for increasing disorder strength (namely for increasing 
a and Uab), while it decreases for increasing interactions 
Ua- This is partly due to the fact that, for a given initial 
distribution (controlled uniquely by the Ua/ Ja ratio and 
by the trapping potential), the gain in on-site interaction 
for the a-bosons is ex Uab — Ua/'i., namely it increases lin- 
early with Uab- On the other hand, increasing Ua reduces 
the probability that many a-bosons occupy the same site 
before the Rabi pulse; this latter fact reduces the prob- 
ability that, immediately after the pulse, a site will be 
occupied by both a and 6-bosons, and this reduces 5e- 

So on the one hand, the equilibrium physics suggests 
that stronger localization (smaller S,) should occur for 
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FIG. 3: (a),(b): transferred energy per a-boson (a) for a non- 
interacting gas of a-bosons (f/a = 0), and (b) for interacting 
a-bosons with Uab = TJa- (c),(d): corresponding localization 
length (in units of the lattice spacing d) of the evolved cloud. 

stronger disorder XJah and for weaker interactions \Ja^ due 
to reduced screening [3|, |7[ • On the other hand, the non- 
equilibrium physics suggests the exact opposite, namely 
that stronger Uah and weaker IJa both increase (5e, so 
that the a-bosons expand more violently, making the 
wings of the quasi-stationary density profile fall off more 
slowly, and therefore increasing the localization length 
^. The evolution of ^ in Fig. [3J:/d shows that the non- 
equilibrium aspects dominate in this competition. Note 
that, for larger ^, the density profile exhibits exponen- 
tial tails over a larger spatial range, which facilitates the 
experimental observation of the tails [18] . 



by disorder more effectively. This effect is reduced upon 
increasing the repulsion strength C/q , showing that disor- 
der screening is indeed observed in the cloud core, which 
is the part contributing more significantly to the overall 
coherence properties. 

In conclusion, we have shown that novel, non- 
equilibrium disorder effects on the many-body state of 
bosons trapped in a spin-dependent optical lattice can be 
observed by dynamically doping frozen particles into the 
system via a Rabi pulse. The doped particles inherit the 
spatial structure and correlations from the initial state 
of the mobile particles, so that a rich disorder statistics 
can be generated by state preparation of the system be- 
fore the pulse. Our predictions can be tested by state- 
selective time-of-flight measurements and in-situ imaging 
of the bosonic cloud, already used to demonstrate An- 
derson localization of weakly interacting Bose gases in 
optically generated random potentials [5|-l7|. Moreover 
an experimental implementation of dynamical doping is 
particularly desirable, as it could access the regime of 
strongly interacting mobile atoms, on the verge of the 
superfluid-to-Mott transition, which remains challenging 
for all current theoretical approaches. 

We thank D. Bauer, J. I. Cirac, and M. Lettner for 
discussions. This work was supported by the European 
Union via SCALA, by the German Excellence Initiative 
via the Nanosystems Initiative Munich, and by the DFG 
via SFB631. 




FIG. 4: Relative coherent fraction R after time evolution (see 
text) for Uab = 7Ja- 



While the real-space localization exhibits counterintu- 
itive, non-equilibrium properties, the coherence proper- 
ties of the expanding cloud exhibit instead a more in- 
tuitive behavior. We consider the momentum distribu- 



tion of the a-bosons, (nfc)t = J2i 



^-ik(t~j)/ J 
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right after the pulse (i = 0) and after time evolution 
to a quasi-stationary state {t = imax = 200?i/Ja), and 
we monitor the loss of coherence during evolution via 
the relative coherent fraction R = {nk=o)t^^^/{nk=o)t=o- 
Fig. H] shows that the relative coherent fraction is sup- 
pressed for increased disorder strength (larger a), show- 
ing that, although the cloud displays exponential tails 
with a longer ^, the internal phase coherence is spoiled 
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